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ABSTRACT 

A  new  technique  for  observer  design  which  requires  the 
solution  of  a  minimax  problem  is  presented.   An  algorithm 
for  solving  the  minimax  problem  is  proposed  and  applied  to 
the  design  of  full-order  and  reduced-order  observers. 
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I.   INTRODUCTION 

In  feedback  control  systems,  the  input  vector  to  a  plant 
is  often  a  function  of  the  current  plant  state  vector.   This 
means  that  one  requires  knowledge  of  the  plant  state  vector 
at  the  current  time  to  accomplish  system  control.   In  most 
systems,  however,  the  entire  state  vector  is  not  available 
for  direct  measurement.   Hence,  it  is  desirable  to  recon- 
struct, or  estimate,  the  unknown  plant  state  vector.   The 
device  which  performs  this  reconstruction  of  the  plant  state 
vector  is  called  an  "observer."   In  this  thesis,  after 
presenting  necessary  background  material  on  the  theory  of 
observers,  which  is  primarily  the  work  of  Luenberger  [4],  a 
new  approach  to  observer  design  is  systematically  developed. 
This  new  method  uses  a  quadratic  performance  measure  which 
is  to  be  minimized  with  respect  to  the  observer  gain  matrix 
G. 

For  simplicity,  only  linear  time-invariant  continuous- 
time  systems  are  considered  throughout  this  thesis. 

The  algorithm  proposed  is  suitable  for  digital  computa- 
tion using  appropriate  subroutines  and  is  applicable  to 
higher-order  systems.   In  order  to  illustrate  the  new 
technique,  numerical  examples  are  presented. 


II.   OBSERVER  THEORY 

A.  THE  NEED  FOR  OBSERVERS 

To  design  a  feedback  control  system  it  is  often  desired 
to  know  the  entire  plant  state  vector.   But,  in  most  cases, 
the  entire  state  vector  cannot  be  measured,  hence  a  suitable 
estimate  to  the  plant  state  vector  must  be  determined.   This 
creates  the  need  for  a  system  called  an  "observer"  which 
produces  an  estimate  of  the  plant  state  vector. 

Assuming  that  it  is  possible  to  build  a  system  which  can 
estimate  some  constant  linear  transformation  of  the  system 
state  vector,  then,  provided  that  the  transformation  is 
invertible,  it  is  possible  to  reconstruct  the  state  vector 
itself.   This  is  the  fundamental  idea  upon  which  observer 
theory  is  based. 

B.  FULL-ORDER  OBSERVERS 

First  consider  the  construction  of  a  full-order  observer 
for  a  free  system- described  by 


x(t)  =  A  x(t)  (1) 

y(t)  «  C  x(t)  (2) 

where  x_(t)  is  the  nxl  state  vector,  y_(t)  is  the  mxl  output 
vector,  A  is  an  nxn  constant  matrix  and  C_  is  an  mxn  constant 
matrix,  and  n  >  m.   Thus,  an  m-dimensional  output  vector  y_(t) 
is  available  from  measurements,  but  we  want  to  know  the 
entire  state  vector  x(t)  of  the  plant.   Since  x(t)  cannot  be 


measured  directly,  one  must  obtain  an  estimate  x(t) ,  that  is 
a  good  approximation  to  the  original  plant  state  x(t)  for 
all  time.   In  what  follows  it  is  assumed  that  the  plant 
described  by  (1)  and  (2)  is  observable. 

One  way  to  estimate  the  state  vector  x(t)  might  be  to 
construct  a  system  of  the  same  dynamic  order  and  structure 
as  the  plant,  given  by 

z(t)  =  A  z(t)  (3) 

and  drive  it  in  parallel  with  the  plant.   If  this  is  done 
the  system  state,  _z(t),  will  be  a  linear  transformation  of 
the  plant  state  x_(t).   If  the  initial  states  z_(0)  and  x_(0) 
are  equal,  then  _z(t)  and  x(t)  remain  equal  for  all  t  £  0. 
Usually,  however,  z_(0)  is  not  equal  to  x(0)  and  the  differ- 
ence between  the  two  states  z_(t)  and  x_(t),  can  be  determined 
by  forming 

z(t)  -  x(t)  =  A  [z(t)  -  x(t)]       (4) 

The  solution  of  (4)  is 

z(t)  -  x(t)  =  e-*  [z(0)  -  £(0)]     (5) 

Thus,  the  error  z_(t)  -  x_(t)  depends  on  the  matrix  A.   If  the 
matrix  A  is  a  stability  matrix,  i.e.,  all  the  eigenvalues  of 
the  matrix  A  have  negative  real  parts,  the  error  approaches 
zero  as  time  goes  to  infinity. 

A  better  alternative  is  to  construct  a  linear  system 
which  contains  a  model  of  the  plant  state  equation  but  is 


also  driven  by  an  output  error  signal,  y_(t)  -  C  z_(t),  that 
is, 

z(t)  =  A  z(t)  +  G[y_(t)  -  C  z(t)]  (6) 


or 

z(t)  =  (A  -  G  C)  z(t)  +  G  y(t)  (7) 

This  system  has  a  very  desirable  feature:  the  observer  gain 
matrix  G  can  be  selected  to  determine  the  rate  at  which  the 
observer  state  z_(t)  approaches  the  plant  state  x(t)  .  If  the 
state  z_(t)  becomes  equal  to  the  plant  state  x(t)  at  some  time 
t  =  t.,  ,  then  the  correction  term,  y_(t)  -  £  ^(t),  is  zero  for 
t  =  t1  and  the  state  z_(t)  will  remain  aligned  with  the  plant 
state  x(t)  for  t  >.  t.  . 

The  system  described  by  equation  (6)  or  (7)  is  referred 
to  as  a  "full-order  observer"  for  the  plant  described  by  (1) 
and  (2).   The  gain  matrix  G  can  be  selected  by  the  designer. 
How  to  select  the. gain  matrix  G  is  the  problem  of  observer 
design. 

The  above  development  can  be  easily  extended  to  forced 
systems  by  including  the  input  u(t)  in  the  observer  equation 
(6)  or  (7).   Consider  the  system  described  by 

x(t)  =  A  x(t)  +  B  u(t)  (8) 

yCt)  -  C  x(t)  (9) 


where  u(t)  is  the  rxl  input  vector  and  B  is  an  nxr  constant 
distribution  matrix.   The  observer  for  this  plant  is  similar 
to  that  for  the  free  system,  but  the  input  vector  u(t)  must 
be  included.   The  appropriate  full-order  observer  system  is 


z(t)  =  A  z(t)  +  B  u(t)  +  G[y(t)  -  C  z(t)] 


(10) 


or 


z(t)  =  (A  -  G  C)  z(t)  +  B  u(t)  +  G  y(t) 


(ID 


The  configuration  of  this  full-order  observer  is  shown  in 
Fig.  1. 


FIGURE  1.   FULL-ORDER  OBSERVER 


C.   REDUCED-ORDER  OBSERVERS 

The  full-order  observer  possesses  a  certain  degree  of 
redundancy.   The  reason  for  this  redundancy  is  that  the  full- 
order  observer  constructs  an  estimate  of  the  entire  plant 
state,  but  the  output  of  the  plant,  which  is  made  up  of 
linear  combinations  of  the  states,  is  available  by  direct 
measurement.   Hence,  it  is  actually  unnecessary  to  construct 
a  system  like  the  full-order  observer  to  estimate  the  entire 
plant  state . 

The  redundancy  can  be  eliminated  by  reducing  the  dynamic 
order  of  the  observer  system  to  (n  -  m) ,  the  difference 
between  the  dimensions  of  the  state  vector  x(t)  and  the 
output  vector  y_(t)  .   The  observer  which  is  reduced  to  its 
minimal  dynamic  order  (n  -  m)  is  called  a  "reduced-order 
observer."   The  structure  of  a  reduced-order  observer  is 
illustrated  in  Fig.  2. 
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FIGURE  2.   STRUCTURE  OF  REDUCED-ORDER  OBSERVER 

The  starting  point  for  consideration  of  reduced-order 
observers  is  again  the  plant  characterized  by 


x(t)  =  A  x(t)  +  B  u(t) 


y_(t)  =  C  x(t) 


(12) 
(13) 
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The  output  distribution  matrix  C_  has  dimension  m.   By  intro- 
ducing a  change  of  coordinates,  it  can  be  assumed  that  the 
matrix  C  takes  the  form  C  =  [  I  J  0  ]  .   The  new  state 
vector  3t(t),  is  selected  so  that  its  first  m  components  are 
equal  to  the  plant  output  vector  y_(t)  and  the  remaining 
(n  -  m)  components,  denoted  by  w(t),  can  be  selected  arbi- 
trarily but  in  such  a  way  that  the  states  y(t)  and  w(t) 
are  linearly  independent.   Thus,  the  new  state  vector  x_(t) 
can  be  written  in  partitioned  form  as 


x(t) 


y(t) 


w(t) 


and  ?[(t)  is  related  to  the  original  plant  state  vector  x(t) 
by  a  nonsingular  matrix  P,  i.e., 


x(t)  =  P  x(t) 


D 


x(t) 


The  matrix  £  here  is  the  one  given  in  (13)  and  D  is  an 
arbitrary  matrix  selected  in  such  a  way  that  P  is  non- 
singular.   The  new  state  equations  have  the  form 


y(t) 
w(t) 


=  A 


y_(t) 
w(t) 


+  B'  u(t)    (14) 


The  matrices  A'  and  B'  in  (14)  are  not  the  same  as  A  and  B 
in  (12).   These  two  sets  of  matrices  are  related  to  one 
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another  through  P_.   Equation  (1*1)  can  be  rewritten  in 
partitioned  matrix  form  as 


y(t) 

A        ' 
-11 

A       ' 
-12 

y_(t) 

■h 

Bi' 

w(t) 

hi ' 

A        ' 
-22 

w(t) 

V 

u(t) (15) 


or 


2_(t)  -  A11t  y(t)  +  A12»  w(t)  +  B1t  u(t) 


(16) 


w(t)  =  A21«  y(t)  +  A22'  w(t)  +  B2«  u(t) 


(17) 


The  vector  v_(t)  is  available  for  measurement,  and  if  it  is 
differentiated,  y(t)  can  also  be  determined.   Since  u(t)  is 
also  available,  (16)  provides  A -,  2 '  w(t)  v;hich  plays  the  role 
of  a  set  of  measurement  equations  for  the  system  (17)  which 
has  w(t)  as  its  state  vector  and  A^  '  v_(t)  +  B_2 '  u(t)  as  its 
input  vector.   The  next  step  is  to  construct  a  full-order 
observer  of  order  (n  -  m)  for  the  plant  of  equation  (17) 
using 

-12*  -(t)  =  ^(t)  "  -ll'  ^(t)  "  -l'  -(t) 

as  the  measurement.  Thus,  the  observer  for  the  system 
described  by  (16)  and  (17)  can  be  expressed  as 
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w(t)  =  (A22!  -  G  Al2')  w(t)  +  A21'y(t)  +  B 2 '  u(t) 


+  1   CZ  Ct)  -  Ai;l'  y_(t)]  -  G  B1'  u(t) 


(18) 


or 


w(t)  =  (A22»  -  G  A12» )  w  +  (A21!  -  G  A11' )  y_(t) 


+  (B2»  -  G  B^')  u(t)  +  G  v_(t) 


(19) 


where  w  represents  the  estimate  of  w.   The  gain  matrix  G 
in  (18)  or  (19)  can  be  selected  so  that  the  observer  system 
coefficient  matrix  App'  -  G  A, p !  has  arbitrary  eigenvalues 
provided  that  (A^p1,  &■??')    is  observable.   It  is  known  [4] 
that  if  (C,  A)  is  observable,  then  (A, 0',  A,*)  Is  also 
observable.   The  configuration  of  the  reduced-order  observer 
is  shown  in  Fig.  3- 
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FIGURE  3-   REDUCED-ORDER  OBSERVER  WITH  DERIVATIVE 
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The  differentiation  of  v_(t)  in  (19)  may  be  troublesome 
to  perform  physically;  however,  it  can  be  avoided  by  modi- 
fying the  observer  state  equation  (19)  .   The  resulting 
modified  observer  state  equation  can  be  written  as 

z(t)  =  (A22«  -  G  A12!)  z(t)  +  (A22«  -  G  A12')  G  y_(t) 


+  (A211  ~  1   Au')  y_(t)  +  (B2'  -  G  B±n    u(t) 


with 


w  =  z(t)  +  G  y(t) 


(20) 


(21) 


and  the  corresponding  observer  diagram  is  shown  in  Fig.  4 
which  is  equivalent  to  the  observer  in  Pig.  3  at  the  point 
w. 


FIGURE  4.   MODIFIED  REDUCED-ORDER  OBSERVER 
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D.   THE  CONVENTIONAL  METHOD  OP  OBSERVER  DESIGN 

It  has  previously  been  shown  that  the  full-order  observer 
for  the  linear  time-invariant  plant 

x(t)  =  A  x(t)  +  B  u(t) 

(22) 
v_(t)  =  C  x(t) 

is  characterized  by 

z(t)  =  (A  -  G  C)  z(t)  +  G  y(t)  +  B  u(t)         (23) 

If  (£,  A_)  is  observable,  the  observer's  eigenvalue  can  be 
located  arbitrarily  by  selecting  the  gain  matrix  G.   The 
problem  is  how  to  select  the  eigenvalue  locations  for  the 
observer. 

The  observer  error  e_(t)  is  defined  as 

e(t)  A  z(t)  -  x(t)  (24) 

Differentiating  equation  (24),  substituting  the  expressions 
for  z_(t)  and  x(t),  and  solving  for  e(t)  gives 

e(t)  =  (A  -  G  C)  e(t)  (25) 

Letting 

P  =  A  -  G  C 

then  the  error  equation  for  the  observer  can  be  written  as 

e(t)  =  F  e(t)  (26) 

15 


The  solution  to  this  equation  is 


e(t)  =  e-  t  e(0)  (27) 


Thus,  the  dynamic  response  of  the  error  between  the  observer 
and  the  plant  state  is  determined  by  the  matrix  F  =  A  -  G  C. 
It  is  desired  that  the  error  e_(t)  eventually  approach  zero 
for  all  e(0).   Thus,  if  G  is  selected  so  that  F  =  A  -  G  C 
is  a  stability  matrix  then  the  observer  state  z_(t)  approaches 
the  plant  state  x_(t)  at  a  rate  controlled  by  G.   One  approach 
to  observer  design  is  to  select  the  desired  eigenvalues  and 
then  determine  G  to  obtain  these  eigenvalues  [*!]. 
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III.   A  NEW  METHOD  FOR  OBSERVER  DESIGN 

A.   PROBLEM  DESCRIPTION 

An  alternative  procedure  for  observer  design  can  be 
based  on  a  performance  measure  for  the  error  system  (26) 
consisting  of  the  integral  of  a  quadratic  form  in  the 
observer  error,  that  is, 

J  =   /    eT(T)  Q  e(x)  di  (28) 

to 

where  Q  is  a  real  symmetric  positive  definite  constant  matrix 
and  e_(i)  is  a  function  of  the  gain  matrix  G,  which  is  to  be 
selected  to  minimize  this  performance  measure.   It  is  shown 
in  the  Appendix  that  equation  (28)  can  also  be  expressed  as 


J  =  e/  V^Q)  e^  (29) 


T 
where  e   is  the  initial  error  of  the  observer,  e   is  the 
— o  — o 

transpose  of  e  ,  and  V  (G)  is  the  value  at  x  =  0  of  the  real 
1        —o}     — o  — 

symmetric  positive  definite  solution  V(t)  of  the  matrix 
differential  equation 

V(t)  =  -FT  V(t)  -  V(t)  F  -  Q  (30) 


subject  to  the  boundary  condition 


V(tf)  =  0  (3D 
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Alternatively,  defining  t  =  tf  -  t,  this  equation  can  be 
expressed  as 


V(t)  =  FT  V(t)  +  V(t)  F  +  Q  (32). 


with  the  initial  condition 

V(0)  =  0  (33) 

From  equation  (29)  it  is  seen  that  J  depends  on  the 

initial  error  of  the  observer,  e  ,  and  on  the  matrix  V  (G) 

'  — o  *  — o  — 

which  is  implicitly  dependent  on  G  through  equation  (30). 
Finding  a  solution  to  equation  (30)  requires  solving 
n(n  +  l)/2  differential  equations  (since  V  is  a  symmetric 
matrix) . 

B.   CONSIDERATION  OF  THE  INITIAL  CONDITION 

It  is  desired  to  select  the  observer  gain  matrix  G  to 
minimize  the  performance  measure 

J  =   /    eT(i)  Q  e(x)  di  =  e  T  V  (G)  e 
—     —  —        — o  — o  —  — o 

To  calculate  the  performance  measure,  a  priori  knowledge  of 
the  initial  state  of  the  error  system  is  needed.   The  initial 
error  state  is  not  known,  but  a  reasonable  approach  is  to 
take  the  worst  initial  state  in  some  bounded  region.   It  will 
be  assumed  that  the  initial  observer  state  is  0_  and  that  the 
possible  initial  error  states  lie  within  a  hyperspherical 
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region  with  center  at  the  origin  and  having  a  specified 
radius  p.   With  this  assumption  the  performance  measure  can 
be  written  as 

fcf    T 

J   =   max    {        f  e  (t)  Q  e(i)  dx 


P 


^o       -   '«■ 


SoNl  P 


— o 


»ax     lej   V^G)  e^  (34) 


llSolli  P 


where  denotes  the  Euclidean  norm.  It  is  easily  shown 
that  the  maximum  occurs  where  |je  ||  -  p  ,  moreover,  without 
loss  of  generality  it  can  be  assumed  that  p  =  1  .  That  is, 
the  maximum  can  be  considered  to  occur  on  the  unit  hyper- 
sphere.  This  is  because,  for  any  e  such  that  lie  1=1  and 
any  radius  p  >  0  , 


J   =  p  e  T  V  (G)  p  e   =  p2  e  T  V  (G)  e 
p  .   — o  — o  —  K  — o    K  — o  — o  —  — o 

=  p2J1  (35) 


i.e.,  the  maximum  performance  measure  on  the  hypersphere  of 
radius  p  is  p  "  times  the  maximum  performance  measure  on  the 
unit  hypersphere.   Hence,  it  is  desired  to  minimize 
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Jl  ■ 


max 


^o 


llsoll 


=  i 


/  e    (t)    Qe(i)    di 


max 


e 

— o 


llfioll 


{ 


=  1 


e    T  V    (G)    e 
— o     — o  — o 


(36) 


with  respect  to  G_. 

C.   EXTENSION  TO  REDUCED-ORDER  OBSERVERS 

So  far,  a  new  method  to  construct  a  full-order  observers 
has  been  discussed.   The  method  can  easily  be  applied  to 
design  reduced-order  observers  also.   Consider  again  the 
linear  time-invariant  system 


x(t)  =  A  x(t)  +  B  u(t) 


(37) 


y(t)  =  C  x(t) 


(38) 


The  reduced-order  observer  equation  for  this  system  is 


z(t)  =  (A22«  -  G  A_12'  )  z(t)  +  (A22»  -  G  A12'  )  G  y_(t) 


+  (A21f  -  1   Au')  y_(t)  +  (B2'  -  G  B^)    u(t) 


(39) 
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with 

w(t)  =  z(t)  +  G  y_(t)  (40) 

where  w(t)  is  the  estimate  of  w(t)  . 

The  characteristics  of  the  observer  system  (39)  depend 
on  the  coefficient  matrix  A_2  '  -  G_  A ,  ?  '  which  depends  on  the 
gain  matrix  G.   Again  to  determine  the  gain  matrix  G, 
consider  the  state  error,  w(t)  -  w(t)  .   The  observer  error 
e_(t)  is 

e(t)  k   w(t)  -  w(t)  (41) 

Differentiating  equation  (^11)  gives 

e(t)  =  w(t)  -  w(t) 
Substituting  w(t)  and  w(t)  and  solving  for  e(t)  yields 

e(t)  =  (A2'2   -G  A12')  e(t)  (42) 

Letting  F  =  A_22'  -  G  A]?',  then  (42)  can  be  written  as 

e(t)  =  F  e(t)  (43) 

The  solution  to  this  equation  is 

e(t)  =  c~  t   e_(Q)  (44) 

Notice  that  equation  (44)  has  exactly  the  same  form  as 
equation  (27)  except  for  the  definition  of  the  matrix  P. 
Therefore,  the  results  obtained  for  full-order  observers 
apply  to  reduced-order  observers  as  well. 
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D.   SIMPLIFICATION  FOR  THE  INFINITE-TIME  CASE 

It  will  now  be  assumed  that  the  upper  limit  on  the 
performance  measure  integral  is  infinity.  It  will  be  shown 
that  if  tf  -*°°  ,  the  problem  of  calculating  V  (G)  can  be 
simplified  to  the  problem  of  solving  a  set  of  n(n  + l)/2 
linear  algebraic  equations.  Let  us  state  an  important 
lemma  and  two  theorems  related  to  this  simplification. 
The  first  step  is  to  write  an  explicit  solution  for  V(t). 

Lemma  1:   The  solution  of  the  matrix  differential 
equation 

V(t)  =  FT  V(t)  +  V(t)  F  +  Q  (45) 


is 

V(t)  =  £(t)  V(0)  £T(t)  +  /  £(t  -  x)  Q  £T  (t  -  x)  dx    (46) 

0 


where  £(t)  is  the  transition  matrix  for  the  system 

z'(t)  =  FT  z(t)  (47) 


Proof:   Suppose  that  (46)  is  the  solution  of  (45),  and  is 
rewritten  as 

V(t)  =  £(t)  V(0)  £T(t)  +  £(t)  [  /  £(-t)  Q  £T(-x)  dx]  ^i) 

0 

k   i(t)  V(0)  iT(t)  +  ±(t)    L(t)  iT(t)  (48) 

where  L(t)  is  defined  as 

L(t)  A  /   £(-x)  Q  iT(-x)  dx  (49) 

0 
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Differentiating  (48)  gives 


V(t)  =  £(t)  V(0)  iT(t)  +  £(t)  V(0)  i(t)  +  i(t)  L(t)  £T(t) 
+  £(t)  L(t)  £(t)  +  £(t)  Ci(-t)  Q  iT(-t)]  £T(t) 


(50) 


The  last  term  of  equation  (50)  becomes  Q,  because 

l(t)  [£(-t)  Q  £T(-t)]  £T(t)  =  £(0)  Q  £T(0)  =  I  Q  I  =  Q 


Also 


i(t)  =  FT  £(t)  (5D 


therefore  (50)  can  be  written  as 


V(t)  =  FT  <J.(t)  V(0)  <j>T(t)  +  <|)(t)  V(0)  <j)T(t)  F 


+  FT  £(t)  L(t)  iT(t)  +  £(t)  L(t)  £T(t)  F  +  Q 


(52) 


But  from  equation  (48) 


£(t)  V(0)  iT(t)  =  V(t)  -  £(t)  L(t)  iT(t)        (53) 

Substituting  (53)  into  (52)  gives 

V(t)  =  FT  [V(t)  -  £(t)  L(t)  £T(t)]  +  [V(t)  -  £(t)  L(t)  £T(t)]  F 

+  FT  £(t)  L(t)  £T(t)  +  £(t)  L(t)  £T(t)  F  +  Q      (54) 
Simplifying  equation  (54)  yields  the  final  desired  result 
V(t)  =  FT  V(t)  +  V  F  +  Q 
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The  following  well-known  theorem  gives  a  necessary  and 
sufficient  condition  for  the  error  system  (26)  to  be  asymp- 
totically stable,  that  is, 


lim  e(t)  =  0        for  all  e(0) 
t->-«> 


Theorem  1:   The  linear  time-invariant  error  system 

e(t)  =  F  e(t)  (55) 

is  asymptotically  stable  if  and  only  if,  given  any  positive 
definite  symmetric  matrix  Q,  there  exists  a  positive 
definite  symmetric  matrix  K  such  that 

FTK+KF=-Q  (56) 

Theorem  1  has  been  proved  in  [2].   This  theorem  can  be 
used  to  determine  whether  or  not  the  error  system  (26)  is 
asymptotically  stable  for  any  choice  of  the  gain  matrix  G. 

The  following  theorem  applies  specifically  to  the  case 
where  tf  ->  °°. 

Theorem  2:   Assume  that  the  linear  time-invariant  system 

e(t)  =  F  e(t)  (57) 

is  asymptotically  stable,  and  has  a  performance  measure 
defined  by 

J  -   /   eT(T)  Q  e(x)  dT  (58) 
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which  is  equivalent  to 


J  =  e  T  V   (G)  e  (59) 

— o  — o   —  — o  \j?j 


where  V(t)  satisfies 


with 


V(t)  =  FT  V(t)  +  V(t)  F  +  Q 


V(0)  =  0 


then 


V(t)  =  FT  V(t)  +  V(t)  F  +  Q  (60) 


has  a  constant  steady-state  solution  V  ,  i.e., 

lim  V(t)  =  V   (a  constant) 
,     —      — o 

t"H» 

Proof:   By  Lemma  1  with  V(0)  =  0_,  the  solution  to  (60)  is 

t  T 

V(t)  =  /   £(t  -  t)  Q  £  (t  -  t)  di  (61) 

0 

T 

TT   +■        m  V  t" 

Substituting  £(t)  -  e—   and  §_   (t)  =  e—   into  (6l)  and 
using  properties  of  the  matrix  exponential,  gives 


rp       4-         m  mm 

V(t)  =  e—    /    e     Q  (e—   e  —   )   dx 
0 

FTt  f      -FTTn   -Ft,    FT  (62) 

=  e—   /   e  —    Q  e  —  dx  e— 

0 
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But  the  system  (57)  is  asymptotically  stable,  so  according 
to  Theorem  1  there  exists  a  positive  definite  symmetric 


matrix  V  such  that 
— o 


Q=-FTV   -V   P 
—     —  — o   — o  — 


(63) 


Substituting  (63)  into  (62)  gives 

V(t)  =  -  e-  t   f      e~-  T  [PT  V  +  V  F]  e"^T  dx  e-* 
0         _  -o   _o  _ 

T 
— "W    t        T 

Since  e  —   and  F   commute,  this  can  be  written  as 


(64) 


T    t  T 

„,,s             Ft  ,  d   r  -Ft  „   -Ft.   Ft 

V(t)  =  -e-   /  --r-Le-  Ve-je- 

0  dT  ~° 


T       T 
Fxt  r  -Ft  „   -Ft 
c—   Le  —   V  e  — 
— o 


1   Ft 
J  e— 


0 


mm  m 

Flt      -Ft  T,  -Ft   Ft     Ft  -  -r   T   Ft 

e—   e     V  e  —  e—  -  e—   I  V   I  e— 

— o  —  — o  — 


T 

=  I  V   I  -  e-  t  V  z-t 
—  — o  —        — o 


(65) 


The  system  (57)  is  stable,  so 


Ft 
lim  e— 


=  0    and 


T 
Ft 

lim   e-  L   =  0 

t-»-°° 


(66) 


Hence,  substituting  (66)  into (65) 


lim   V(t)  =  V    (a  constant) 

,     —      — o 

t-x» 


which  proves  the  theorem. 
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According  to  the  result  of  Theorem  2,  V(t)  becomes  zero 

as  t  -*■ »  .   Therefore,  V  is  the  solution  of 

— o 

PTV+VF=-Q  (67) 


which  is  obtained  by  setting  V(t)  =  0_  in  equation  (32). 

Equation  (67)  gives  a  set  of  n(n+l)/2  linear  algebraic 

equations  that  can  be  solved  for  V  after  any  positive  definite 

Q  has  been  selected  for  the  performance  measure  (28).   Hence, 

the  problem  of  calculating  V  becomes  simple.   The  positive 

— o 

definiteness  of  the  solution  V  can  be  checked  by  using 
Sylvester's  criterion  [7]  and  if  and  only  if  V  is  positive 
definite  the  system  is  asymptotically  stable. 
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IV.   SOLUTION  OF  THE  OBSERVER  DESIGN  PROBLEM 

As  mentioned  before,  the  error  equation  needed  to 
determine  the  gain  matrix  G  has  the  same  form  for  both  the 
full-order  and  reduced-order  observers.   So  one  design 
method  can  be  applied  for  both  cases. 

A.   DEVELOPMENT  OF  THE  METHOD 

It  is  desired  that  the  gain  matrix  G_  be  found  which 
minimizes 

T 
J_  =   max    \      S     e    (x)  Q  e(x)  dx 


1 


0 


e 

— o 


e   =  1 
— o  " 


max     [  e  T  V  (G)  e   ]  (68) 

— o  — o  —  — o 

e 

— o 

ell  -  1 


Letting  min  J,    =   J-.*>    then 


J   *   =  min          max  [e   T  V    (G)  e    ]                        (69) 

1                  r,  — O        — O  — O 

2-  e 

— o 


llsoll-  i 
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This  is  a  minimax  problem,  and  equation  (69)  is  a  real 

quadratic  form  in  e  .   There  is  a  useful  theorem  about  the 
^  — o 

extremal  properties  of  the  eigenvalues  of  a  real  quadratic 
form  [1], 

Theorem  3:   The  global  maximum  of  a  real  quadratic  form 
on  the  unit  hypersphere  is  equal  to  the  largest  eigenvalue 
of  the  quadratic  form,  and  moreover  the  corresponding  eigen- 
vector is  the  vector  drawn  from  the  origin  to  the  point  on 
the  hypersphere  where  the  quadratic  form  achieves  its 
maximum. 


Therefore 


max      [  e  T  V  (G)  en  ]  =  A  (G)  (70) 

— o  — o  —  — o      1  — 

e 
— o 

•oil"  X 


where  X, (G)  is  the  largest  eigenvalue  of  V  (G).   Hence,  now 
1  —  o  — 

the  problem  is  to  minimize  X-.CG)  with  respect  to  the  elements 
of  G;  that  is,  to • find  the  observer  gain  matrix  G  such  that 


J  *  =  min   X  (G)  (7D 

l     Q    1 


To  ensure  that  the  observer  is  stable  and  does  not  have  an 
unlimited  bandwidth,  the  gain  matrix  G  will  be  selected  from 
the  set  JJ   defined  as 

1.   All  eigenvalues  of  F  have  negative  real  parts.' 


r         1.   All  eigen- 
2.   Sum  of  eij 


;envalues  of  F  =  trace  F  >  -D 
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where  D  is  a  specified  positive  constant.   For  a  given  value 
of  G_  the  largest  eigenvalue  is  easily  found  by  using  a 
digital  computer.   The  program  required  is  simple  because 
very  efficient  methods  are  available  for  finding  the  largest 
eigenvalue  of  a  real  symmetric  matrix,  and  many  computer 
facilities  have  several  standard  subroutines  that  can  be 
used  for  this  purpose. 

B.   ALGORITHM  FOR  COMPUTING  THE  GAIN  MATRIX 

The  steps  in  the  computational  procedure  for  finding  the 
optimal  gain  matrix  are  given  below. 

1.  Guess  G  subject  to  the  following  requirements: 

a.  F  is  stable.   This  can  be  done  by  assigning  arbitrary 
eigenvalues  with  negative  real  parts  to  F  and  using  Luenberger's 
method  [4],  for  example. 

b.  Trace  F  _>  -D 

2.  Use  a  minimization  subroutine  to  minimize 


J..  =  max      ]  e  T  V  (G)  e  \ 
1  .  \   — o  — o  —  — o 

IISoll-  1 


,)-., 


Actually,  a  penalty  function  is  added  to  ensure  that  the 
trace  F  _>  -D  constraint  is  satisfied,  e.g., 

r  0,  if  trace  F  _>  -D 

Jo   =  Jl        +        \  P 

a    x         (.a*  [trace  F  +  D]  ,  if  trace  F  <  -D 

=  Ji  +  JP 

where  a  is  a  positive  weighting  factor. 
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In  computing  values  of  J   for  the  minimization  subroutine 

d 

the  following  steps  are  carried  out. 

a.  Solve  F  V  +  V  F  =  -  Q  for  V  =  V 

b.  Check  to  see  that  V  is  positive  definite.   If  it 

— o 

is  continue,  otherwise  set  J   =  10  ,  that  is,  a  large 
positive  number  and  return  to  the  minimization  subroutine 

c.  If  V^  is  positive  definite  it  is  the  steady-state 
solution  of 

V(t)  =  FT  V(t)  +  V(t)  F  +  Q 


Next  find  the  largest  eigenvalue  of  V  using  a 
suitable  subroutine. 

d.  Calculate  the  trace  of  F  and  compute  the  penalty 

function  J  . 
P 

e.  Set  J   =  J.,  +  J  and  return  to  the  minimization 

alp 

subroutine . 


In  applying  this  algorithm  to  the  examples  given  in  this 
thesis  the  minimizations  were  accomplished  using  a  subroutine 
"DIRECT"  (NFS  Computer  Facility)  that  performs  the  pattern 
search  method  of  Hooke  and  Jeeves,  and  the  largest  eigen- 
values were  found  using  a  subroutine  "GIVHO"  (NPS  Computer 
Facility)  that  utilizes  the  Givens-Householder  method. 
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V.   APPLICATIONS  OF  THE  ALGORITHM 

A.   A  FULL-ORDER  OBSERVER 

The  purpose  of  this  example  is  to  illustrate  the  appli- 
cation of  the  described  algorithm  to  the  design  of  a  full- 
order  observer  for  a  second-order  system. 


Y(S) 


Consider  the  plant  vT ,  Q  \    =  —~ with  the  state- 

iHbj    s^  +  3S  +  2 


variable  representation 


x(t)  =  A  x(t)  +  B  u(t) 


0    1 
■2   -3 


x1(t) 
x2(t) 


u(t) 


(72) 


y(t)  =  C  x(t) 


[1-   0] 


x1(t) 

x2(t) 


(73) 


A  full-order  observer  for  this  system  is  determined  by 
specifying  the  observer  gain  matrix 


g- 


G  = 


32 


The  observer  system  coefficient  matrix  is 


P  =  A  -  G  C  = 


-«1    1 
-2-g0   -3 


(7*0 


which  has  the  characteristic  equation 


+    (3  +  g±)X   +    3&±   +   g2   + 


2  =  0 


(75) 


Selecting  the  observer  eigenvalues  arbitrarily  at  —3  ±  j  3 
gives  the  characteristic  equation 


( A  +  3  +  j  3)  (X  +  3-J3)  =  X  +  6x  +  18  =  0 


(76) 


Equating  the  coefficients  of  like  powers  of  X  in  (75)  and 
(76)  yields 

3 


G  = 


7 


This  is  the  starting  guessed  value  of  the  gain  matrix  G.   If 
one  arbitrarily  chooses  Q  =  I_3  and  solves  (67)  for  V  the 
result  is 


V 


15  +  Sg-j^  +  5g9  +  &2' 


11 


V 


12 


V 


22 


12  +  22gl  +  6g2  +  2glg2  +  6&1' 

3  -  2Sl  -  glg2 ^ 

12  +  22g;L  +  6g2  +  2gxg2  +  6gl' 

3  +  3gl  +  g2  +  gx2 ; 

12  +  22gl  +  6g2  +  2g;Lg2  +  6E;L' 


(77) 


33 


Due  to  the  symmetry  of  V  only  three  equations  are  needed. 
If  the  sum  of  the  observer  eigenvalues  is  to  be  limited  to 
-10.0,  for  example,  then  the  proposed  algorithm  gives  the 
solution  for  the  observer  gain  matrix  as 


G  = 


7.0 
■0.64063 


(78) 


which  gives  the  observer  eigenvalues  as 


X  - 


-6.625 
-3.375 


(79) 


For  other  selections  of  the  weighting  matrix  Q  and  different 
limiting  values  for  the  trace  of  F  the  values  obtained  for 
the  observer  gain  matrix  G  and  the  corresponding  eigenvalues 
are  tabulated  in  Table  1.   Looking  at  the  table  it  can  be 
seen  that  the  two  eigenvalues  become  farther  apart  from  each 
other  when  the  limiting  value  of  the  sum  of  the  eigenvalues 
is  decreased  (assuming  a  fixed  Q) . 

To  compare  this  full-order  observer,  which  will  be 
referred  to  as  the  optimal  observer,  with  other  observers 
having  different  observer  eigenvalues,  consider  two  observers 
having  eigenvalues  selected  as 


X  = 


-5.0 
-5.0 


X  = 


-8.0 
-2.0 


(80) 


respectively.   The  corresponding  observer  gain  matrices  are 
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G 


7.0 
2.0 


G   - 


7.0 
•7.0 


(81) 


respectively. 

For  simulation  purposes,  the  input  u(t)  was  selected  to 
be  a  unit  step  function.   Assuming  that  the  initial  error 
states  lie  on  the  unit  circle,  the  actual  and  the  estimated 
states  computed  by  using  the  optimal  gains  of  equation  (78) 
are  plotted  as  functions  of  time  on  the  graphs  in  Fig.  5  for 
the  best  initial  error  condition  and  in  Fig.  6  for  the  worst 
initial  error  condition.   The  estimate  errors  for  the 
observers  with  three  different  observer  eigenvalues  are 
plotted  on  one  graph  in  Figs.  7  and  8  for  the  best  and  worst 
initial  error  conditions,  respectively.   The  values  of  the 
performance  measure  for  the  observer  error  systems  with 
different  observer  eigenvalues  (different  observer  gains) 
are  tabulated  in  Table  2  for  the  worst  initial  conditions. 


Observer  Gains 
(Eigenvalues) 


Worst  Initial 
Conditions 


Values  of  Perform- 
ance Measure 


Si  =  7-0 


Dl 


% 


Al  = 


6.62b' 


=  -0.64063Un  =  -3.375, 


e(0)  = 


0.16 
-0.987 


0.164 


&i   = 


7.0 
2.0 


'K   =  - 


A^  =  - 


5.0 

5.0 


e(0)  - 


0.525 
-0.85 


0.181 


&i 


S2 


7.0 
-7.0 


A,  =  - 


A„  p  - 


8.0 
2.0 


e(0)  = 


0.63  ' 
0.777 


0.303 


TABLE  2.   PERFORMANCE  VALUES  FOR  VARIOUS  OBSERVER 
EIGENVALUES 
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B.   A  REDUCED-ORDER  OBSERVER 

The  purpose  of  this  example  is  to  illustrate  the  appli- 
cation of  the  algorithm  to  the  design  of  a  reduced-order 
observer  for  a  third-order  plant. 

Consider  the  plant  characterized  by  the  transfer  function 


Y(S)  _ 


U(S)    S3  +  4S2  +  5S  +  2 


This  plant  can  be  represented  by  the  state  and  output 
equations 


xx(t) 

x2(t) 
x3(t) 


0    10 

0    0    1 

•  2   -5   -4 


x1(t) 
x2(t) 
x-(t) 


+ 


u(t) 


(82) 


y(t) 


[1   0   0] 


xx(t) 
x2(t) 
x3(t) 


(83) 


This  plant  is  a  third-order  system  with  a  single  output  so 
a  second-order  observer  with  arbitrary  eigenvalues  can  be 
constructed.   The  output  matrix  C  already  has  the  same  form 
as  the  one  developed  in  the  reduced-order  observer  theory. 
Hence,  it  is  not  necessary  to  change  the  variables,  that  is, 
for  this  case  the  matrices 


A'  =  A 


B'  =  B 
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The  matrices  A  and  B  can  be  partitioned  as 


A 


0  ,'  1 

0  |  0 
» 

-2  «  -5 


0 

1 
-4 


B  = 


Therefore,  the  reduced-order  observer  system  matrix  is 


I  =   A22  -  G  A12 


&1 
-5-g, 


1 


(84) 


which  has  the  corresponding  characteristic  equation 


X*  +  (4  +  g1)X  +  4g1  +  g2  +  5  =  0 


(85) 


Suppose  one  selects  the  starting  observer  eigenvalues  at 
X  =  -3j  -2.   This  gives  the  characteristic  equation 


(X  +  3)(X  +  2)  =  X  +  5X  +  6  =  0 


(86) 


Equating  the  coefficients  of  the  same  powers  of  X  in  (85) 
and  (86)  yields 


G  = 


1.0 
•3.0 


This  is  the  initial  guessed  value  of  the  observer  gain  matrix 
G.  If  the  weighting  matrix  Q  is  chosen  arbitrarily  as  Q  =  I, 
and  equation  (67)  is  solved  for  V  the  result  is 
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V 


36  +  4g;L  +  llg2  +  g^ 
40  +  42gl  +  8g2  +  2glg2  +  8gl2 


V 


4  "  5gl  "  glS2 


40  +  42gx  +  8g2  +  2glg2  +  8g;L2 


(87) 


6  +  4gl  +  g2  +  gl: 


22   40  +  42g;L  +  8g2  +  2glg2  +  8gl2 


Due  to  the  symmetry  of  V  only  three  equations  are  needed. 
The  proposed  algorithm  gives  the  solution  for  the  observer 
gain  matrix  as 


G  = 


26.0 
•1.73438 


(88) 


if  the  sum  of  the  eigenvalues  of  F  is  limited  to  -30.0. 
This  gain  matrix  G_  gives  the  observer  eigenvalues 


A   = 


-25.851 
-4.149 


(89) 


To  compare  this  optimally  designed  reduced-order 
observer  with  others  having  different  observer  eigenvalues, 
consider  two  observers  with  the  eigenvalues 


4  3 


X   = 


■15.0 
•15.0 


-20.0 
-10.0 


(90) 


respectively.   The  corresponding  observer  gain  matrices  are 


G 


26.0 
116.0 


G   = 


26.0 
91.0 


(91) 


respectively.   In  Fig.  9  the  actual  and  estimated  states 
using  the  optimal  gains  are  plotted  for  the  best  initial 
error  conditions  and  in  Fig.  10  for  the  worst  initial  error 
conditions.   It  is  assumed  that  the  initial  error  condition 
lies  on  the  unit  circle.   The  estimate  errors  for  the  three 
observers  are  plotted  on  one  graph  in  Figs.  11  and  12  for 
the  best  and  worst  initial  error  conditions  respectively. 
Table  3  shows  the  values  of  the  performance  measure  for 
observers  with  different  eigenvalues  for  the  worst  initial 
conditions. 

C.   OBSERVATIONS 

The  examples  indicate  that  the  observers  designed  by 
using  the  proposed  technique  avoid  large  errors  at  a  cost  of 
having  relatively  small  errors  that  slowly  approach  zero. 
This  is  to  be  anticipated  from  the  form  of  the  performance 
measure  selected. 

The  design  approach  is,  however,  conservative  in  the 
sense  that  the  resulting  observers  function  effectively  for 
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Observer   Gains 
(Eigenvalues) 


Worst   Initial 
Conditions 


Values   of  Perform- 
ance  Measure 


gx  =  26.0 


'A     =  -25.851^ 


^  =  -1.731138  \\2  ■-     -4.1119, 


e(0) 


0.122 
-0.993 


0.123 


&i    = 


% 


26.0 
116.0 


fx1  =  -15.0 

A2  =  -15.0 


e(0)  = 


0.98 
-0.2 


1.150 


g2 


26.0 
91.0 


'^1  =  -20.0 
^  =  -10.0 


e(0)  = 


0.966 
-0.26 


0.842 


TABLE    3.       PERFORMANCE   VALUES    FOR   DIFFERENT   OBSERVER   EIGENVALUES 
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situations  in  v/hich  the  initial  conditions  are  v/orst.   In 
exchange  for  this  characteristic,  performance  for  less 
difficult  initial  conditions  is  compromised. 
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V.   CONCLUSIONS 

In  this  thesis  a  new  method  for  designing  observers  has 
been  developed.   The  optimal  observer  gains  are  determined 
by  solving  a  minimax  problem.   The  proposed  algorithm  is 
easily  suited  to  digital  computation  and  especially  appli- 
cable to  higher-order  systems. 

The  observers  computed  by  using  the  algorithm  have  a 
tendency  to  avoid  large  errors  at  the  expense  of  small 
errors  which  slowly  decay.   As  anticipated,  the  observers 
designed  by  the  proposed  technique  are  efficient  for  the 
worst  initial  error  conditions. 
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•  ~  ~  I"'1 ' 


•:j;ere    e  [z 


_r.=    s  :  _  . 


:: 


, 


c  -        -  —    -  — 


-      N      —     */■*- 


-    z)   e(t)  (A.«) 


=  : 


and  to  the  initial  state  by 

e(tf)  =  £(tf)  e(0)  (A. 5) 

Differentiating  (A. 4)  with  respect  to  t  gives 

i(tf  -  t)  =  -  £(tf  -  t)  F  (A. 6) 

for  all  t~  and  t  with  the  boundary  condition 
&(tf  -  tf)  =  £(0)  =  I 

The  performance  measure  for  (A.l)  is  given  by 

tf 
J  -  /    [  eT(x)  Q  e(x)  ]  dx  (A. 7) 

0 

Substituting  for  e_(i)  using  (A.  2)  gives 

J  -  /     [  eT(0)  £T(t)  Q  1(t)  e(0)  ]  dx  (A. 8) 

0 


Rearranging  (A. 8),  one  obtains 

fcf 

J  =  eT(0)  [  /   £T(t)  Q  £(t)  dx  ]  e(0)  (A. 9) 

0 


Defining 

tf 

V(t)  £  /    £T(t  -  t)  Q  £(t  -  t)  dx  (A. 10) 

t 

then        . 
f 
V(0)  =  /    ^T(t)  Q  1(t)  dx  (A. 11) 

0 
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Therefore,  we  have  the  final  result  (29) 


J  =  eT(0)  V(0)  e(0)  (29) 


Differentiating    (A. 10)    with   respect    to   t   by   using  Leibnitz's 
rule 

fcf 
V(t)    =    /  [   4)T(t   -   t)    Q  £(t   -   t)    +  £T(x   -    t)    Q  £(x   -   t)    ]    dx 


-   |T(t    -    t)    Q  i(t   -    t) 

f  [   £T(t   -   t)    Q  ±(t   -   t)    +  £T(x   -   t)    Q  £(x   -    t)    ]    dx 

t 


-   Q  (A. 12) 

Substituting  (A. 6)  for  _£(  x  -  t)  gives 
t 


V(t)  =  /    {  -  [PT  £T(x  -  t)  Q  £(t  -  t)] 
t    i 

-  [^T(t  -  t)  Q  £(x  -  t)  F]|   dx   -   Q 


(A. 13) 


Rearranging  (A. 13)  gives 

V(t)  =  -  FT  [   /    ^T(t  -  t)  Q  £(x  -  t)  dx] 

t 

fcf 

-  [   /    4T(t  -  t)  Q  £(x  -  t)  dx]  F  -  Q        (A. I'D 
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Substituting  for  V(t)  using  (A. 10)  gives  equation  (30) 


V(t)  =  -  FT  V(t)  -  V(t)  P  -  Q  (30) 


Letting  t  =  tf  in  (A. 10)  gives  equation  (31) 


V(tf)  =  0  (3D 


Taking  the  transpose  of  (A. 10),  the  right  side  does  net 
change,  hence  V(t)  is  symmetric. 


55 


BIBLIOGRAPHY 


1.  Bellman,  R. ,  Introduction  to  Matrix  Analysis,  McGraw-Hill, 
I960. 

2.  Kalman,  R.  E.  and  Bertram,  J.  E.,  "Control  System 
Analysis  and  Design  via  the  'Second  Method'  of  Lyapunov: 
I.   Continuous-Time  System,"  Journal  of  Basic  Engineering, 
Trans.  ASME,  Series  D,  Vol.  82,  pp.  371-393,  I960. 

3.  Kirk,  D.  E.,  Optimal  Control  Theory:   An  Introduction, 
Prentice-Hall,  1970. 

k.      Luenberger,  D.  G.   "An  Introduction  to  Observers,"  IEEE 
Trans.  Automatic  !  ntrol,  Vo .  AC-16,  pp.  596-602,  Dec. 
1971. 

5.  Ozer,  E.,  "Design  of  Constrained  Optimal  Controls  for 
Linear  Regulators  with  Incomplete  State  Feedback,"  E.E. 
Thesis,  N.P.S.,  Monterey,  California,  1969. 

6.  Perkins,  W.  R. ,  and  Cruz,  J.  B. ,  Engineering  of  Dynamic 
Systems,  John  Wiley  &  Sons,  Inc.,  New  York:  1 969. 

7.  VJiberg,  D.  M.,  Schaum' s  Outline  Series:  State  Space  and 
Linear  Systems,  New  York:   McGraw-Hill,  1971. 

8.  Wilf,  H.  S.,  Mathematics  for  the  Physical  Sciences,  John 
Wiley  &  Sons,  Inc.,  New  York:   1962. 


56 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia   2231^ 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California  939^0 

3.  Professor  Donald  E.  Kirk,  Code  52Ki  5 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   9  39^0 

4.  Professor  Sydney  R.  Parker,  Chairman  1 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   939^0 

5.  LCDR  Chae-Young  Park,  Korean  Navy  5 
Naval  Headquarters 

Young  Dung  Po,  Seoul 
KOREA 


57 


SECURITY   CLASSIFICATION   OF   THIS  PAGE  (When  Dili  Entered; 


REPORT  DOCUMENTATION  PAGE 

HEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.     REPORT  NUMBER 

2.  GOVT  ACCESSION  NO. 

3.     RECIPIENT'S  CATALOG  NUMBER 

4.     TITLE  (and  Subtitle) 

Time-Domain   Design   of  Observers 

5.     TYPE  OF   REPORT  &  PERIOD  COVERED 

Master's    Thesis; 
December   1973 

6.     PERFORMING  ORG.   REPORT  NUMBER 

7.     AUTHORf*; 

Chae-Young  Park 

C     CONTRACT  OR  GRANT  NUMBERfij 

9.     PERFORMING  ORGANIZATION   NAME   AND  ADDRESS 

Naval   Postgraduate   School 
Monterey,    California      939^0 

10.     PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

It.     CONTROLLING  OFFICE  NAME   AND   ADDRESS 

Naval   Postgraduate   School 
Monterey,    California      939^0 

12.     REPORT  DATE 

December   1973 

13.     NUMBER  OF   PAGES 

59 

14.     MONITORING  AGENCY  NAME  &    ADDRESSf*/  dlffvtcnt  from  Controlling  Oillce) 

Naval   Postgraduate   School 
Monterey,    California      939^0 

IS.     SECURITY  CLASS,  (of  this  report) 

Unclassified 

ISe.    declassification/downgrading 

SCHEDULE 

16.     DISTRIBUTION   STATEMENT  (of  this  Report) 

Approved    for  public   release;    distribution  unlimited. 

17.     DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  II  different  from  Report) 

IB.     SUPPLEMENTARY  NOTES 

19.     KEY  WORDS  (Continue  on  reverse  aide  It  neceeeiuy  and  Identify  by  block  number) 

Pull-Order  Observer 
Reduced-Order   Observer 
Observer   Gain 
Observer  Error 
Minimax 

20.     ABSTRACT  (Continue  on  reverse  tide  If  necessary  and  Identify  by  block  number) 

A   new   technique    for   observer   design   which   requires    the 
solution   of   a  minimax  problem   is   presented.      An  algorithm   for 
solving   the   minimax  problem   is   proposed   and   applied   to   the 
design   of   full-order   and   reduced-order   observers. 

0D     1  JAN   73     1473  EDITION  OF    1   NOV  65  IS  OBSOLETE 

(Page      1)  S/N    0102-011-6601   | 

58 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Bntured) 


4'UCUHITY   CLASSIFICATION   OF   THIS  PAGEfWhan  Date  Entered) 


DD     Form       1473  (BACK) 

1  Jan  73  , 

S/N      0102-014-6601  t-g        SECURITY  CLASSIFICATION  OF  THIS  PAGEf»?ion  Due  Ent»r*d) 


P1573  Park 
c.l 


Jfme-domafn  desi 
observers. 


gn  of 


Thesis 
P1573 

c.l 


Park 

Time-domain   design  of 
observers. 


0*6^6 


thesP1573 

Time-domain  design  of  observers. 


3  2768  001  00047  4 

DUDLEY  KNOX  LIBRARY 


